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Abstract 



In a recent paper, we showed that tree level up-down unification of fermion 
Yukawa couplings is a natural consequence of a large class of supersymmetric 
models. They can lead to viable quark masses and mixings for moderately 
large values of tan /3 with interesting and testable predictions for CP violation 
in the hadronic sector. In this letter, we extend our discussion to the leptonic 
sector focusing on one particular class of these models, the supersymmetric 
left-right model with the seesaw mechanism for neutrino masses. We show 
that fitting the solar and the atmospheric neutrino data considerably restricts 
the Major ana- Yukawa couplings of the leptons in this model and leads to 
predictions for the decay r fi -\- j, which is found to be accessible to 
the next generation of rare decay searches. We also show that the resulting 
parameter space of the model is consistent with the requirements of generating 
adequate baryon asymmetry through lepton-number violating decays of the 
right-handed neutrino. 



The discovery of neutrino mass has added a new impetus to the search for new physics not 
only beyond the standard model but also beyond the minimal supersymmetric extension of 
the standard model (MSSM). One candidate for such a theory is the supersymmetric version 
imQ] of the left-right symmetric model , which provides a natural embedding of the seesaw 
mechanism for small neutrino masses. There are two ways of realizing the seesaw mechanism 
in these models: One by using the B — L = 2 iso-triplet Higgs fields in which case the seesaw 
arises from purely renormalizable interactions; another using B — L = 1 iso-doublet Higgses 
to get seesaw out of nonrenormalizable interactions. The first alternative has the merit 
that it results in automatic i?-parity as a consequence of gauge invariance. The second 
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alternative has the apparent advantage that known compactification schemes of heterotic 
string models can lead to this class of models at low energies. The disadvantage is that 
one needs extra symmetries to maintain the known degree of baryon and lepton number 
conservation (equivalently i?-parity). 

In a recent paper we have shown that in both these versions of the SUSY left-right 
models (as well as in their generalizations to 5*0(10)), if the theory below the right handed 
SU{2)ji scale is assumed to be given by the MSSM, then one ends up with a strong tree 
level relation between the up and down quark and lepton Yukawa coupling matrices, i.e., 

Yu =Y,; Y, = Y,n. (1) 

We call this up-down unification 0. These relations follow for arbitrary Higgs structure 
in the first class of models (i.e. the ones with B — L = 2 HIggs fields) whereas in the 
second class of models they hold only in the minimal version with one bi-doublet field. 
Furthermore, in the second class of models, there is an overall constant of proportionality 
instead of equality. In any case, as is clear, these relations lead to vanishing of the quark 
mixing angles (the CKM angles) at the tree level. They also lead to the proportionality 
predictions for the quark masses (m^ : : = : rris : rrih) which do not agree with 
observations. As was demonstrated in Ref. [^], once the one loop corrections coming from 
gluino exchanges (which are known to be significant in such models, especially when tan/3 
is large [^]) are included, the observed mixing and masses for the quarks can be obtained for 
a particular choice of flavor mixing in the squark sector. It was shown in Ref. that the 
flavor changing neutral currents that are associated with such squark mixings are consistent 
with experimental limits. Further, it was pointed out that these models lead to specific 
predictions for the CP- violating parameters e'/e and electric dipole moment of the neutron 
which are different from the conventional MSSM so that they can be used to test these 
models. Additional tests of these models are provided by branching ratios of Higgs decay 
into fermion pairs, and deviation from the standard model predictions in rare processes such 
as Bfi — Bd mixing. 

It is the goal of this paper to explore the implications of up-down unification in the lepton 
sector. The interest in this problem stems from the fact that up-down unification reduces the 
number of parameters in the seesaw formula for the neutrino masses without the assumption 
of grand unification. For instance, in a generic seesaw model for three generations, there are 
twelve (complex) parameters describing the complete seesaw matrix for the light neutrinos 
(three from the heavy Majorana neutrino sector and nine from the Dirac mass matrix). If 
one assumes that the neutrino Dirac mass is symmetric (as can happen in some GUT models 
or left-right models), we are still left with 9 parameters. On the other hand, in models with 
up-down unification that we discuss in this paper there are only six parameters. This is 
because the Dirac neutrino masses in a particular basis is diagonal and completely given by 
the charged lepton masses. This narrows down the possible parameter space in which a joint 
fit to the solar and atmospheric neutrino data can be obtained. We find that the Majorana 
Yukawa couplings that fit data exhibit a hierarchical pattern similar to that in the charged 
fermion sector of the standard model. Such a hierarchical structure in the heavy Majorana 
neutrino mass matrix goes well with the idea of baryogenesis through lepton asymmetry 
induced by the decay of the lightest right-handed neutrino {Nf.). We use the neutrino fits 
to make predictions for the rare decay r — /i + 7 as a function of the seesaw scale (denoted 
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by vn) . Such rare decays have a supersymmetric origin: In running the soft sfermion masses 
from the scale of gravity (Mpianck or Mstrmg) to vr, the right-handed neutrinos are "active" 
through their Majorana Yukawa couplings, which induce flavor mixing in the slepton sector. 
We find that there exists a range of vr below the GUT scale for which the predicted r — /i+7 
branching ratio is within the reach of the planned experiments. 

Let us start by giving a brief derivation of the up-down unification relation in the SUSY 
left-right models. As is well known, the electroweak gauge group for this model is SU{2)l x 
SU (2)/j X U (1) B-L with the standard assignment where Q, denote left-handed and right- 
handed quark doublets, denote the (2,2,0) Higgs bi-doublets, A and A*^ are the left- and 
right-handed Higgs triplets and A and A^ fields, their conjugates. Let us write down the 
gauge invariant part of the matter superpotential involving these fields: 

+ i{iL^T2AL + feL'^^TaA^L'^). (2) 

Left -right symmetry implies h^*| = (h^*|)^ and f = fp. It is clear from the above equation that 
if there is only a single bidoublet field $(2,2,0), then up-down unification will follow. As 
has been shown in Ref. [0], this result holds even when more than one bi-doublet is included 
in the theory. The result is more general if i? — L = 2 Higgs fields are employed for SU{2)r 
gauge symmetry breaking. The same result holds if the right handed symmetry is broken 
hy B — L = 1 doublets as well, but if there is only one Higgs bi-doublet. What happens in 
this case is that at the vr scale, the two doublets (0^,^) from the bi-doublet $ and the two 
left-handed doublet superfields transforming as (2, 1, ±1) under SU (2) l x SU {2)r x U (1) b-l 
(call them x and x) form a 2 x 2 mass matrix. Once we assume that appropriate fine tuning 
has been done to leave the MSSM doublets (denoted by Hu and Hd) light below the vr 
scale, their generic form becomes = cos 6'0„ + sin 6xu and Hd = cos a(j)d + sin axd- In this 
case, we obtain the modified up-down unification formula F„ = 7!^ and = 'jY^d, where 
7 = cos 6/ cos a is a proportionality constant. 

Let us now study the implications of up-down unification in the lepton sector. We pass 
to a basis where the charged lepton mass matrix is diagonal. In this basis one gets: 

M^D = tan p Mi (3) 

for the renormalizable seesaw case (triplet case) and the same relation with an extra constant 
multiplier in the doublet case. Here, Mi = Diag{mf., m^, 711^-). The light Majorana neutrino 
mass matrix is then given by: 

tan^ 3 

Vr 

where f is the right-handed neutrino Majorana Yukawa coupling defined in Eq. (1). In the 
case of using doublet Higgs fields, f happens to be the coupling associated with the non- 
renormalizable interaction of the right handed neutrinos i.e. fij{Lfx'^ L'jx'^) / M and analogous 
discussions will follow. Eqs. (2)- (4) parameterize the neutrino masses and mixings in terms 
of six parameters. All parameters are complex in general, but for simplicity of presentation 
we shall assume them to be real. 
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To study the neutrino masses and mixings in this model, let us start with a specific 
parameterization of the / matrix as follows: 




(5) 



This is the most general form of a symmetric 3x3 matrix, except for the in the (1, 3) 
entry. This choice has been motivated by considerations of baryogenesis, discussed later in 
this article, which requires /13 to be very small. Eq. (5) leads, after see-saw diagonalization, 
to a Majorana neutrino mass matrix given by: 

= / r X i,N jmerrif, -m^ m^mM ■ (6) 

vn9{e^ 6, A, b) ^ _^.^^2 _ ^) ) 

Here g{e,5, A, 6) = y — 6^(1 — 6^). Note that the neutrino masses and mixings depend on 
four independent parameters which are completely fixed by observations, thereby fixing all 
parameters of the lepton sector of the theory. Also note the suppression factors arising from 
the charged lepton masses in Eq. (6). 

Recently, the Super-Kamiokande collaboration ^ has established a deficit in the mea- 
sured flux of muon neutrinos from the atmosphere. This observation can be interpreted as 
compelling evidence in favor of neutrino masses and oscillations. The most likely scenario 
is oscillation of —>■ 1/^. The mass-splitting is inferred to be: Am^^ ^ (10^^ — 10"^) 
eV^ and the oscillation angle is in the range: sin'^26^r — 0.8 — 1. We will assume param- 
eters consistent with this interpretation in our analysis. The long-standing solar neutrino 
deficit can be explained by neutrino oscillation in matter (MSW oscillation between 
z/g and another species (most likely z/^ if there are no extra species of neutrinos, as in our 
models). Two possible solutions exist: the small angle MSW and large angle MSW re- 
gions. We shall adopt the small angle solution for which Am^^ ^ (0.3 — 1) x 10^^ eV^ and 
2 X 10~^ < sin^ < 2 x 10~^. To see how one obtains solutions that fit both solar and 
atmospheric neutrino data, let us first focus on the — sector and assume that ^ ^. 
In this case, it is clear that if e ~ — , then for 6 ~ 1 we have maximal mixing between the 
z/^ and i>r. Furthermore, if 6 = 1, we have rriij = 0. Therefore if we take b = 1 + n with 
K ^ 1, we get 

^ - - (7) 

m,^ 2 ^ ^ 

and the z/^ — z/^ mass difference sqared given by: 

9 9 tan'^ 3 (4:b)mt 



We thus see that mj^^ <^ m^^ and the atmospheric neutrino data fixes m^^. To get solar 
neutrino via the MSW small angle solution, we require k ^ 0.1 or so. The formula for m^^ 
then enables us to determine the value of vpt and we find it to be around 10^^ GeV if we 
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choose tan (3 ~ 40, which seems required from considerations of the quark mixings. Let us 
hasten to note that this is not a complete determination of the seesaw scale - since we have 
scaled the Majorana coupling f^^ = 1 to get this value. There is an overall scaling freedom 
and we could scale all the fi/s by a common factor that would change the v^. For instance, 
if we scaled down all the /j/s by a factor of 100, we would get vr ~ 10^'' GeV which is near 
the GUT scale. If we limit the value of vn to stay below the string scale (~ 4 x 10^'' GeV), 
then we have a small range for /ss between 1-0.01. Thus, although there is no dynamical 
mechanism in the model that would independently fix the overall scale of the /j/s, we do 
not have too much freedom in their range. 

The solar neutrino puzzle has also a solution in this specific parameterization. To see 
this, note that one easily finds niu^ <^ m^^ and we have adjusted the value of n to get the 
right mass difference squared between i/^ and u^. The mixing has a suppression factor of 
rrie/m^ but is multiplied by the ratio 6/\ so that we can adjust those parameters to get the 
required mixing angle needed in the small angle MSW solution. We do not mean to imply 
that this is a unique solution - but rather we present this example to illustrate the general 
nature of the detailed numerical fits using the parameterization for the f matrix. Below we 
give the values for the different elements of the f matrix that lead to acceptable solutions 
to both the solar and the atmospheric neutrino data. 

With tan /? = 41 and vr = 10^^'^ GeV we choose 

7 X 10-^ 

-2.5 X 10-3 4.8 X 10"=^ I ; (9) 




4.8 X 10 



~2 




which gives rise to neutrino masses: (4.3 x 10 ^, 2.9 x 10 ^, 6.1 x 10 ^ ) eV and the 
leptonic mixing matrix: 

-0.99 

U = \ 2.8 X 10-2 0.64 -0.77 I ; (10) 




-3.8 X 10-2 

In our notation, U12 is the mixing angle relevant for solar neutrinos. 

We have also obtained consistent fits including non-zero /13. For example, choosing 

/-l.lxlO-5 -1.0x10-^ 3.8xl0-3\ 
f= -1.0x10-^ -1.7x10-3 4.8 X 10-2 . (11) 
V 3.8 X 10-3 4.8 X 10-2 -1.56 / 

gives rise to following neutrino masses at v^i = 10^^'^ GeV (for tan/3 = 41): (9.8 x 10-^, 3.0 x 
10~3^ 6.4 X 10"2) eV. The leptonic mixing matrix is given by: 



U = 



/ -0.99 

9.4 X 10-2 -0.70 0.71 I . (12) 

V3.2 X 10 



-2 




There are small corrections to these values coming from the renormalization group extrap- 
olation from the vr scale to the weak scale ||11[| . 



Let us now comment on the baryogenesis constraint on the matrix elements of f. In 
the picture where baryogenesis results from the leptogenesis caused by the decay of heavy 
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Majorana neutrinos in the very hot stage of the universe |T^, one needs these neutrinos to 



be out of equihbrium. It is well known that in these scenarios, only the lightest of the three 
right handed neutrinos (the A^e) is responsible for leptogenesis. In the presence of /13, its 
decay has a contribution coming from the r + final state and is given by 

. (13) 

Furthermore, the mass of the N^. must be less than the re-heating temperature of the universe 
after inflation, which is required by gravitino abundance constraints to be less than 10^ GeV. 
Putting this value in the out of equilibrium condition r(A^'e) < H requires that /13 < 10^^. 
Thus the /13 element must indeed be very small and this justifies our choice of zero for this 
parameter. With the form of the heavy Majorana neutrino matrix fixed from low energy 
neutrino oscillation data, we can examine if the required lepton asymmetry can be induced 
in the decay of Ne [0. We find that adequate lepton asymmetry, as large as 10~^, can be 



easily induced with the choice of parameters given in Eqs. (9), (11). The actual baryon 
asymmetry will be quite sensitive to the dilution factor. Given that a lepton asymmetry of 
order 10~^ can be realized in our model, some dilution arising from N^. not totally being out 
of equilibrium will be allowed. This is the reason that in the second example (Eq. (11)) we 
allowed a small /13 entry which is marginally inconsistent with the bound in Eq. (13). 
Predictions for Rare r and ^ decays 

Having obtained the various elements of the matrix f , let us turn to the phenomenological 
implication of this result in the leptonic sector. We shall assume that the slepton masses 
are universal at the Planck or GUT scale. Now due to the presence of off-diagonal elements 
in f , the renormalization group evolution will induce mixings between the second and the 
third generation sleptons (i.e. terms such as /i^f ) as well as between the first and the second 
generation sleptons. Such mixings arise because between the Planck/GUT scale (where 
universality holds) and vr, the Majorana Yukawa couplings will contribute to the evolution 
of the soft parameters. The jl — f mixing will depend on the product /33/23. Similarly, there 
will be mixing between e and fi proportional to /11/12 as well as /13/23. These off diagonal 
slepton mass terms are well known to induce rare lepton decays such as r ^ /i + 7 and 
/i ^ e + 7. 

There are two possible diagrams in lepton flavor violating decay modes in this model: (i) 
the chargino (x^) mediated (involving neutral sleptons) and (ii) the neutralino {x^) medi- 
ated (involving charged sleptons). The chargino diagram is comparable with the neutralino 
diagram. In the chargino diagram, we have contributions only from the left-handed sneutri- 
nos. This is because only the left-handed sneutrinos survive down to the weak scale. On the 
other hand the left as well the right-handed charged sleptons contribute to the neutralino 
diagram. We will see that the masses of both the helicities develop flavor violating soft 
terms through RGE due to the presence of off diagonal elements in f . 

Ij —^li + '~f has the following effective lagrangian {mi. > mi.): 

L,ff = !!^{Aj,a^''PRl,F,, + A^ka^'^PJ^F,,) + h.c. (14) 

where Ai = Al^± + Al^o and A2 = A2^± + A2^o. 
The decay width is given by: 
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nh^ki) = ^^{\M' + \M\ (15) 

Al^± and Al^± are given by: 

^V^ m=lfc=l 



'J 



= -wj^, E E X {i/r/'^i^:r'='^i^2(x,.,j} , (17) 



where the convention Xab = ""^a/^^b has been adopted and G^j^'-' = fC*™'^r^^ if^''^'* = 
y^(jrn,i ^^j^ifj^ ' . are the sneutrino rotation matrices and the uC, vC are the chargino 

rotation matrices and Ui = diag{me,m^,mr) /(V^MwcosP). M^± are the chargino masses 
and m?^ are the sneutrino masses. -^2(3^) = l/(12((a; — l)'^)(2x^ + Sx"^ — 6a; + 1 — Gx'^lnx) 
and Fsix) = l/i2{{x - lf){x^ - 4a; + 3 + 21nx). 

The expressions for the neutrahno contributions are [O]. 



-^^--w^££sr^ '''' 



{(v/2GS?) (v^GSf)f.ks,;)- 

771/ ~ 

/ \ / iilj 



. ay/a ^ ^ 1 



rrij 



where 



G-f = -l/2[Z„i + cot ^H^Z^^lrfi: (19) 

Gmkh y -pkh 

OIR — ^mli iJi 

kh 



Tjmkh r7 / -p j-T 

rlniT — Zjm-A i U.'J] 



OIL - {i^lL^l 

-mki ry /-p JT 

OIR — '^rnS U IRUi 



where Z is the 4x4 neutrahno mixing matrix in the (^B, W3, H^, ifg) basis, F4^{x) = 1/ (2((x— 
l)^)(x^ — 1 — 2xlnx) and Il,r. = ^Ilr^- slepton mass eigenstates. /^^/j is 6 x 6 
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slepton mass matrix l^^n at the weak scale and can be written in the 3x3 forms having the 
submatrices Mll,Mlr and Mrr as follows: 

/ MilMI Mf^^ \ 



where 



{MilMI),, = (ral^).. + M| cos 2(5 Q - sin^ Ow) Sij , (21) 



Here Ui is the lepton Yukawa matrix, Ui^ = AiUi and the A's are the trilinear couplings 
present in the potential. We determine the soft masses and the trilinear couplings at the 
weak scale by using the RGE, which will be discussed next. 

Before we calculate the branching ratios of r — > /i + 7, let us first discuss the parameter 
space of this model. We assume the model originates at the string scale (Mstring ~ 10^^ 
GeV) and is left-right symmetric. The trilinear terms involving the A fields, bi-doublets 
and the matter fields in the superpotential are shown in Eq. (2). We observe that the A 
and A'^ have couplings given by f and fc to the leptons (fc = f due to left-right symmetry). 
This superpotential is valid upto the scale vr. At this scale all the A fields and the right- 
handed neutrinos pick up mass. The doubly charged components of A and Ac fields do 
not acquire masses of order vr, their masses arise through non-renormalizable operators 
and are of order ^Ij/Mpianck- This situation is different if we employ SU{2) singlet doubly 
charged Higgs fields, which may be desirable from the point of view of inducing electroweak 
symmetry breaking (see comments later). The RGE of the soft supersymmetry breaking 
slepton masses are given by: 

^ = ^[-47r(3/2«B-LM^_^ + 3«^Ml) (22) 

+ \{{UiU} + f ft)mL + raldUiU} + f ft) + 2{Uml^U}) 

+ 2(fmi^ft + mlUiU} + m^f ft + UiJJI + f//,f/)j)]. 

The same RGE also holds for the right-handed sleptons because of left -right symmetry. It 
is easy to see that the terms like f ftm|^, Uf^Uj^ (f//^ = Aft) introduce fiavor violations 
in the soft terms for the left and the right handed sleptons. However, the amount of fiavor 
violations that is introduced through Uf^Uj^ are rather small. 

The RGE's for the other soft terms, Yukawa couphngs and the trihnear terms are as 
follows: 

dm^ 2 

^ = 16^ [-4vr(3ai? A4 + 3aiM|) (23) 
+ 3A?(2mJi + mJ + A|) 
+ 2Tr{mlUiUl + Uinil^Ul)], 
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^ = Y^[-4^(4«L^1 + 4«b-lM|_J (24) 
+ 2Tr{mlUiU} + t/,mLt// + C//,t//A^)], 

^ = + ^"^^ + ^«s-i)t^iA (25) 

+ 87r(3a^Affl + MlMl + ^as-LMB_L)[/L 

+ 8C/,^ + AUiUlUi^ + 4f ftf/,^ + 2UfJjM 

+ 2Tt{UiJJI)Ui + Tr{UiUl)Ui^ + SAtX^Ui + SXpui 

dUf 19 9 

= I^t-^^C^"^ + 2'^B-L)Uf^ + 8TT{7aRMR + -aB-LMB-L)f (26) 

+ 8f/,,f//f + AUiUjUf, + AUfjH + 2f ftf/^, 
+ 2Tr(^/,ft)f + Tr(fft)t/^J, 

f = T^[-^^(^"« + \ocb-l)1 + 2f ft + 4C/,C// + Tr(f ft)f], (27) 

^ = l^t"^''^^''^ + ^""^ + + 2f ft + 4t/;.t// + Tr{3UgU,^ + t/^[/,t)t/;]. (28) 

Uq is the up quark (or the down quark) Yukawa matrix. Because of the left-right symmetry 
an, q;l, f and renormahze identically. Mis arc the gaugino masses. 

The doubly charged component of the A"^ fields, which got leaked through, picks up 
mass at the scale Ma(~ jt^ — )• The superpotential, between Vji and Ma, contains a new 

-''^string 

interaction term (in addition to the MSSM interactions) given by i^l'^A'^^^l'^, where /'^'s are 
the charged leptons. The right-handed slepton masses are affected by this interaction term 
and the RGE is given by: 

^ = ^[-4.(12M) (29) 

+ iml + f^^)inl^ + ml^ml + f^t) + 2{Uiml^Ul) 
+ {rmlj^r^ + 2ml^UiU} + mLiif'^ f'^^ + 2UiJJI + UfiU}.J)]. 

There is a further twist to the story of soft masses in our model. The constraint At = 
at the left-right scale renormalizes the up type and the down type Higgs identically. It is 
almost impossible to make the (mass)^ of one of the Higgses to be positive and the other 
negative at the weak scale, unless there exists some splitting of the up and down type Higgs 
masses already at the left-right scale. Since in this model both the up type and down type 
Higgs are coming from the same bi-doublet, there is no scope to introduce non-universality 
in their masses at the string scale. However, there exists some natural splitting for these 
masses at the left-right scale arising from the D term. These D terms get generated due to 



9 



the breaking of SU{2)r x U{1)b-l U{1)y and depend on the mass difference of the A'^ 
and the A'^ fields at the left-right scale. For example, the up type Higgs mass-squared gets a 

correction: —-rrr—^ — — ^ and the down type Higgs gets a negative of the above correction. 

If the masses of A'^ and A'^ are the same at the string scale, m\c will be smaller than m^^ 
due to renormalization, since Ac has Majorana Yukawa couplings. The mass-squared of the 
down type Higgs will then get a negative contribution while the up type Higgs will get a 
positive contribution. This goes against electroweak symmetry breaking requirement. The 
sign of this D term splitting is fixed if universal boundary condition is assumed at the string 
scale. We can however assume that the soft breaking masses for the A"^ and the A'^ fields 
are not same (at the string scale) and that these masses are are also different from the other 
soft breaking masses. A rationale for this assumption is discussed later. We denote the 
departure from universality in the masses by 5j. We assume that the soft mass squared term 
for the A"^ field at the string scale is given by + Si) and for the A"^ field is given by 
ml{l + 62) ■ TTLo is the soft breaking mass of the other scalars. With this assumption, we find 
that the electroweak symmetry can be broken in a large region of parameter space when 
we assume Si to be positive and S2 to be negative. This non universality in the A masses 
will also percolate down to other scalar masses at the vn scale through the D terms. The 
lighter slepton masses gets lowered and this increases the BR(r fi + '-f). We will assume 
the gaugino masses to be universal. 

One possible source for the non-universality in the masses of A'^ and A at vr is that it 
might be dynamically generated by Yukawa couplings. As noted earlier, with the minimal 
set of Higgs fields, the doubly charged components of A'^ and A'^ acquire masses of order 
w^/Mpianck- One way to give these members masses of order vr is to introduce (1,1, ±4) 
Higgs fields {rj^rf). Then the following Yukawa couplings in the superpotential are allowed: 
AiA^A^r] + A2A A 77. If A2 < Ai, at vr, the soft mass-squared of A will be smaller than 
that of A'^. The D term splitting will then have the right sign to lower the mass of Hy^ 
relative to that of if^, facilitating electroweak symmetry breaking. 

For most part of our numerical analysis, we shall assume a spectrum without the 77,7/ 
fields. We have verified that the major consequence of including these fields can be taken 
care of by allowing non-universal soft masses for the A'^ fields at the Planck scale. One 
example for this case is given in Fig. 6. 

Using these masses as input at the string scale, we determine the spectrum of this model 
at the weak scale. The procedure of computation goes as follows. We first determine the 
Yukawa couplings and f at the string scale using the experimental masses and mixing angles 
of the quarks and leptons. At the string scale we introduce the soft breaking masses and 
run down all the masses and the couplings to the left-right scale using the RGBs relevant to 
the left-right symmetry. From the scale, at which the doubly charged fields get decoupled. 



down to the weak scale, we use the MSSM RGEs(as given in ref. (we are also using 
their sign convention of the /i term). In between the left-right and the doubly charged field 
decoupling scale we use the MSSM RGEs coupled with the Eq. (|30|) . At the weak scale we 
have the 6x6 slepton mass matrices Eq. ([20|) . All the 3x3 blocks are evaluated at the weak 
scale for the charged slepton masses. For the neutral sleptons, M;^^ and the Mj^MiR get 
decoupled at the left-right scale. 

We also determine the other sparticle masses which are not much different from what 
could be their masses in MSSM for a given value of mo, mi/2, tan/? and A . 
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The constraint = A;, at the left-right scale has further restrictions. The equality 
of the Yukawa coupling makes it impossible to have an experimentally allowed mass 
[nib = 4.25 ± 0.2) for low and intermediate values of tan f3. We find that tan /5 ~ 36 — 45 is 
allowed by the experimental constraint for reasonable values of S. The one loop correction 
is also included for the evaluation of the rrib and this loop contribution involves a gluino 
mediated diagram (involving sbottom) and the chargino mediated diagram (involving stops). 
The one loop correction is given by [|T^,0: 

Arrib ~ -[^^M^fitanPg{ml^,ml^, M^) + j^^^Atfitanf] g^fht^.mt^, i/)]\bVd (30) 

where g{a, 6, c) = J+^^^'|^°^Ll+° a ] ) ^^^^ jg ^j^g gluino mass, mfe. are the sbottom 

and rht^ are the stop masses. The sign of /i is needed to be positive (in our convention) in 
order to have correct value of irib. But it has been also realized pr5|-[T7| that for this choice 



of the sign of /x, satisfying the CLEO bound (1 x 10"^ < BR(6 ^ s + 7) < 4.2 x 10"^ pl) 
on BR(6 — s> s + 7) requires larger values of mo and mi/2. The reason for this is that the 
chargino diagram adds constructively to the SM and the Charged Higgs mediated diagrams. 
In Fig |1], we plot the BR(6 — >■ s + 7) as a function of mi/2- We have shown curves for 
different values of itiq and 5,. We find that there exists parameter space for mg > 600 and 
mi/2 > 600 which is allowed by this rare decay. 

In Fig ^ we plot the total mass of the 6-quark (mf,), which includes the one-loop cor- 
rection, as a function of mo. We have drawn curves for two different values of mi/2 for 
tan/? =41. 

We now present our calculations of BR(r ^ /i + 7) in the allowed regions of parameter 
space. The experimental upper limit on the BR(r ^ /i + 7) mode is 3 x 10~^ at the 



90 % C.L. [T^. This is expected to improve in the near future [Q. In Fig. y, we plot 
k = Logio ^^^IqIg^'^^ as a function of mi/2- We show curves for different values of mo. We 
find that, as 61 increases and 62 decreases, the BR(r — >• + 7) gets enhanced (the reason 
is that the righthanded and the left-handed stau mass gets lowered). We also find from 
the figures that the BR(r — > yU + 7) is within one or two order of magnitude below the 
experimental value. In this model we also have — > e + 7 mode, but the branching ratio is 
too small (almost 5-6 order of magnitude below the experimental observation). The reason 
for this is that the first and second generation Majorana Yukawa couplings are much smaller 
than the third generation ones. Consequently, the first and second generation slepton masses 
are not as suppressed as the third generation slepton masses. 

So far we have given our results for universal trilinear term A = 0. If we assume A to be 
non-zero, we can lower the rrib mass correction (since the chargino contribution increases the 
mass correction). The allowed parameter space of this model would increase. However one 
should be careful about the lighter stau mass. The change in At at the weak scale is much 
less than the change in Ar for a non-zero Aq. Consequently the lighter stau mass-squared 
(which is already light in this model) can turn negative. However one also can choose A 
judiciously, so that the stau mass is lighter (but not lighter than the lightest neutralino) and 
as a result the BR(r — yU + 7) will increase. We show examples of such case in Figs ^ |^. 
In these figures we can see that the BR(6 — s -|- 7) does not change at all compared to the 
Fig ^ but the BR(r — > /i + 7) changes quite a bit. 
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In the Fig. ^ we plot the BR(r — > /i + 7) in the case when the doubly charged Higgs fields 
get decoupled at the left-right scale. This is the situation that would arise in the presence 
of rj, rj fields, as discussed earlier. The slepton mass matrix will have less fiavor violation in 
this case and the Br{T — /i + 7) is smaller than the previous cases. The Br{h ^ s + 7) 
remains unchanged since the squark sector does not couple to the doubly charged fields. 

In conclusion, we have shown that in the left-right models for neutrino masses, the 
up-down unification allows a fit to both the solar and atmospheric neutrino data with a 
hierarchical set of Majorana Yukawa couplings of the leptons. In the resulting parameter 
space, one obtains a prediction for the rare decay r — >^ + 7 which is within the accessible 
range of the planned experiments. 

The work of KSB is supported by funds from the Oklahoma State University. RNM is 
supported by the National Science Foundation grant No. PHY-9802551. 
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FIG. 1. BR of 6 ^ s + 7 is plotted as a function of the universal gaugino mass (?Tix/2) for 
tan/? = 41. The dotted line is for Si =3, S2 =-0.9 and tuq = 600 GeV. The dashed line is for 
5i =2.5, 62 =-0.9 and mo = 700 GeV. The dash-dotted hue is for 61 =2, 62 =-0.9 and tuq = 800 
GeV. 
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FIG. 2. m5(m{,)is plotted as a function of the universal gaugino mass (^1/2) for tan/3 = 41. 
The parameters for the curves are same as in fig. 1 
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FIG. 3. BR of r — > // + 7 is plotted as a function of the universal gaugino mass (mi/2) for 
tan/3 = 41. The parameters for the curves are same as in fig. 1 
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FIG. 4. BR of 6 — > s + 7 is plotted as a function of the universal gaugino mass (w-1/2) for 
tan/? = 40. The dotted hue is for =3, 82 =-0.9, mo = 650 GeV and ^0 = -550 GeV. The 
dashed line is for 5i =3, 82 =-0.9, mo = 700 GeV and Aq = —600 GeV. The dash-dotted line is for 
5i =2.5, 62 =-0.9, mo = 800 GeV and Aq = -650 GeV. 
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FIG. 5. BR of r ^ ^ + 7 is plotted as a function of the universal gaugino mass {mi/2)- The 
parameters for the curves are same as in fig. 4 
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FIG. 6. BR of T — > /i + 7 is plotted as a function of the universal gaugino mass (^1/2) for 
tan /? = 41 in the case when doubly charged fields get decoupled at the left right scale. The dotted 
line is for 5i =3, 62 =-0.9 and mo = 600 GeV. The dashed line is for 61 =2.5, 62 =-0.9 and ttiq = 700 
GeV. The dash-dotted hue is for 61 =2, 62 =-0.9 and mo = 800 GeV. 
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